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BRIEF:  Image  depth  profiling  is  used  to  determine  relative 
sputtering  rates  of  features  in  plant  tissue  analyzed  by 
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Ion  microscopy,  a  form  of  secondary  ion  mass  spectrometry  (SIMS) , 
is  an  analytical  technique  with  the  ability  to  generate  mass  resolved 
images  of  a  solid  sample  with  elemental  detection  limits  at  the  parts 
per  million  level  (1) .  A  2-20  keV  primary  ion  beam  is  used  to  sputter 
the  top  atomic  layers  of  the  sample  and  secondary  ions  emitted  from 
this  region  sure  extracted  and  mass-analyzed  to  produce  ion  images. 
Qualitative  interpretations  of  ion  images  have  been  shown  to  be  useful 
in  studies  of  elemental  distributions  in  biological  tissues  (2-8) .  However 
quantitative  interpretations  of  ion  images  can  be  erroneous  because  of 
the  problem  of  differential  sputtering,  a  phenomenon  first  pointed  out 
by  Galle  (9)  and  more  recently  by  Parmer,  et  al  (10) . 

The  observed  SIM S  secondary  ion  intensity  may  be  described  by  the 
equation  of  Morrison  and  Slodzian  (1) : 

I  -  TC  Si  a  (1) 

m  p  o 

*<«*....«  I  is  the  measured  secondary  ion  current,  t  the  practical  ion  yield, 

C  the  atomic  concentration  of  the  element  M  (corrected  for  its  isotopic 

A 

abundance),  S  the  total  sputtering  yield  (number: of  atoms  of  any' Kind 

removed  per  incoming  particle) ,  i  the  number  of  primary  particles 

P 

arriving  per  unit  time  and  unit  area  on  the  target  surface,  and  aQ  is  the 
analyzed  area  on  the  target  surface. 

Local  variations  in  composition  across  a  spatially  heterogeneous 
sample's  surface  cause  different  regions  of  the  sample  to  sputter  at 
different  rates.  Since  the  observed  intensity  I  is  proportional  to  S, 
this  causes  erroneous  ion  intensity  enhancements,  and  therefore  higher 
apparent  elemental  concentrations,  for  the  more  rapidly  sputtered  regions. 
Biological  tissue  is  a  very  heterogeneous  material,  and  is  thus  particu¬ 
larly  prone  to  differential  sputtering.  . 


An  appropriate  correction  for  differential  sputtering  is  therefore  needed 
for  any  type  of  quantitative  work  for  biologica)  tissue  in  particular,  and 
spatially  heterogeneous  materials  In  general. 

This  paper  demonstrates  a  first-order  differential  sputtering  correction 
using  digital  image  processing.  It  is  performed  by  mounting  a  thin  biol¬ 
ogical  sample  upon  a  substrate  which  produces  secondary  ions  at  a  mass  not 
observed  in  the  sample.  Images  at  masses  of  interest  are  recorded  from 
the  sample.  The  sample  is  then  sputtered  down  to  the  substrate  while  a 
series  of  ion  images  is  recorded,  with  increasing  depth,  at  the  substrate’s 
mass  by  the  image  depth  profiling  (IDP)  method  of  Patkin  and  Morrison  (1i). 

A  "burn-through11  map  is  generated  from  the  IDP  images,  recording  the  time 
each  location  of  the  sample  sputtered  down  to  the  substrate.  This  time- 
domain  burn-through  map  is  then  used  to  correct  the  initial  sample  ion 
images  for  differential  sputtering  induced  ion  intensity  artifacts. 


EXPERIMENTAL  SECTION 


Instrumentation:  Ion  images  were  obtained  with  a  CAMECA  IMS-3f  ion 
microscope.  This  ion  microscope  obtains  images  at  selected  mass/charge 
ratios  which  retain  the  original  spatial  relationships  of  the  elements  in 
the  sample- (12) . 

The  microscopic  image  digital  acquisition  system  (MIDAS)  ( 1 3)  shown  in 
Fig.  I  was  used  to  record  images.  MIDAS  consists  of  a  low-light  level  I S IT 
(intensified  silicon  intensified  target)  T.V.  camera,  digital  frame  buffer 
(digitally  stores  a  256  x  240  x  12-bit  image),  anaiog/dig itai  converter,  computer, 
graphics  display  screen,  and  associated  computer  software.  The  T.V.  camera  is 
focused  upon  the  fluorescent  screen 
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of  the  ion  microscope.  The  images  are  acquired  by  the  T.V.  camera/frame 
buffer  and  stored  by  the  computer.  Changes  from  the  original  MIDAS 
description  (13)  include  replacement  of  the  CAMECA  IMS-300  ion  microscope 
with  more  sophisticated  CAMECA  IMS-3f  and  dedicated  H.P.  9845T  micro¬ 
computer,  replacement  of  the  DEC  PDP-11/20  minicomputer  and  GT-40  graphics 
display  terminal  with  m  DEC  PDP-11/34A  with  256K  byte  RAM  and  two  10M  byte 
disks,  and  the  use  of  a  f/2.0,  102-mm  zoom  lens  on  the  T.V.  camera  set 
at  40  mm  for  a  SO  yn  field-of-view. 

Computer  Software:  Programs  were  written  in  FORTRAN  IV,  RATFOR 
(a  structured  FORTRAN  IV  pre-processor) ,  and  MACRO- 11  assembly  language. 

There  are  four  stages  for  differential  sputtering  correction  by 
image  depth  profiling  and  digital  image  processing;  acquisition  of  the 
image  data,  amendment  of  the  data  for  detector  and  SIMS  induced  artifacts, 
image  display,  and  the  calculation  of  the  relative  differential  sputtering 
rates  of  various  sample  features.  Five  separate  but  related  programs 
were  used  for  these  functions. 

IMAGE  performs  T.V.  image  acquisition,  disk  input/output  of  images, 
image  display  and  registration. 

GRFTDK  acquires  a  series  of  images  and  transfers  them  to  disk 
storage  at  an  operator  defined  rate  of  up  to  one  image  every  eight  seconds. 

DETCOR  corrects  images  for  spatial  inhomogeneities  in  the  image 
detector  system.  This  is  performed  by  acquiring  a  standard  image  of  a  smooth 
Si  <100>  single  crystal,  with  the  IMS-3f  immersion  lens  optics  slightly 
defocused  to  provide  a  uniformly  intense  ion  image.  The  MIDAS  output  image 
is  influenced  by  non-uniform  detector  response  and  degradation  so  that  it  appears 

inhomogeneous.  DETCOR  determines  the  correction  factor  at  each  pixel 
necessary  to  bring  that  pixel's  intensity  up  to  the  maximum  intensity  of 
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the  standard  image,  correcting  for  detector  system  heterogeneity.  It 
then  uses  this  factor  to  correct  the  current  elemental  image. 


GCONV  converts  raw  image  intensities  to  the  corresponding  ion 
intensity  using  an  empirically  derived  calibration  equation. 

DISCOR  generates  the  bum-through  map  and  then  uses  it  to  correct 
the  image  for  differential  sputtering.  This  is  done  by  initially 
setting  each  pixel  of  the  map  to  zero.  The  substrate  images  acquired  by 
GRFTDK  are  then  sequentially  analyzed  and  the  bum-through  map  updated 
as  follows: 
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where  1^  ^  is  the  new  value  of  the  pixel  located  at  x,y  in  the  bum-through 
map,  current  value  is  the  value  currently  stored  in  that  location  of  the 
bum-through  map  and  next  value  is  the  intensity  at  location  x,y  of  the 
image  currently  being  examined  in  the  GRFTDK  depth  profile  sequence.  T  is 
the  sputtering  time  at  which  the  GRFTDK  image  being  examined  was  acquired. 

Simply,  this  functions  as  follows:  Since  the  substrate  on  which 


the  sample  is: mounted  is  being  imaged,  any  pixel  with  a  non-zero  intensity 
indicates  the  sample  has  been  completely  sputtered  away  at  that  location. 
To  avoid  spurious  readings,  however,  the  program  waits  until  a  location 
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gives  a  continuous  series  of  non-zero  substrate  ion  intensities  in  the  IDP 


sequence  before  assigning  that  location  a  bum-through  time.  The  bum- 
through  time  assigned  is  the  time  at  which  the  first  of  these  images  was 
acquired.  This  procedure  is  repeated  for  every  image  in  the  IDP  sequence, 
ultimately  yielding  the  bum-through  map  which  contains  the  time  at  which 
each  pixel  first  sputtered  through  to  the  substrate. 

At  this  point  the  bum-through  map  is  used  to  perform  a  linear 
correction  to  the  observed  image  intensities  of  the  elements  of  interest 
in  the  sample  such  that: 


„c  i  bum-through  time 

I  =■  I  .  - 3 - 

X'^  X total  sputtering  time 

c  i 

where  I  is  the  differential  sputtering  corrected  image,  and  I  is  the 

detector  discrimination  corrected  and  ion  intensity  converted  elemental  image 

This  gives  an  ion  image  corrected,  at  least  to  a  first  approximation, 

for  differential  sputtering  effects.  This  correction  is  based  on  two 

assumptions.  The  first  is  that  the  observed  secondary  ion  intensity  is 

proportional  to  the  sputter  rate.  The  second  assumption  requires  that  a 

given  region  of  the  sample  sputters  at  a  uniform  rate  throughout  the  course 

of  the  analysis. 

Analytical  Procedure:  The  IMS-3f  instrumental  conditions  are  shown 
in  Table  I.  The  parameters  were  chosen  to  optimize  image  quality,  intensity, 
and  sample  sputtering  rate.  The  time  between  images  was  chosen  on  the  basis 
of  the  primary  ion  beam  species  and  sample  thickness. 
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Sample  Preparation:  Root  tfps  of  Raphanus  satfvus  (radish) 
seedlings  were  selected  for  the  study  because  of  their  welt  defined 
intracellular  calcium  distribution.  One  to  2-mm  pieces  of  root  tips 
were  chemically  fixed  for  1  hr.  using  5%  gluteraldehyde  in  0.03  M 
Pipes  buffer  [piperazine  N-N^b is, (2-ethanol  sulfonic  acid) ]  at  pH  6.8  (14). 
Root  tips  were  washed  with  0.15  Pipes  buffer  and  then  postfixed  for  1  hr. 
in  1%  osmium  tetroxide  in  the  same  buffer.  After  dehydration  in  ethanol 
(70-90$)  and  propylene  oxide,  the  tissue  pieces  were  embedded  in  low 
viscosity  resin  (15).  Sections  were  cut  on  distilled  water  using  a  LKB 
Ultrotome  III  and  transferred  onto  a  polished  tantalum  disk.  Tantalum  was 
selected  as  a  substrate  because  it  was  not  expected  to  be  present  in 
tissue.  Sections  were  then  stored  under  vacuum  until  analyzed  in  the  ion 
microscope. 

RESULTS  AND  DISCUSSION 

Calcium  in  plant  root  tip  cells  was  chosen  for  analysis  due  to  its 
well  defined  intracellular  distribution  (16).  A  light  micrograph  of  the 
stele  cell  region  of  a  radish  root  tip  section  is  shown  in  Figure  2. 

Figure  3  is  a  40Ca+  ion  micrograph  clearly 

showing  higher  intensity  of  calcium  in  nuclei  and  cell  walls  than  in 
cytoplasm.  The  severity  of  differential  sputtering  is  indicated  by 
Figure  4.  This  is  not  an  image  of  40Ca+,  but  of  the  substrate,  lslTa+, 
taken  in  the  course  of  the  image  depth  profile.  These  features  clearly 
show  that  differential  sputtering  causes  the  nuclei  and  cell  walls  to  sputter 
away  well  before  the  cytoplasm  region.  This  is  a  serious  effect  which 
must  be  considered  in  any  quantitative  work. 
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Figure  5  is  a  typical  burn-through  map,  with  the  lighter  areas  indicating 
the  locations  requiring  the  longest  time  to  sputter  away.  The  map  clearly  shows 
the  cell  walls  as  the  most  rapid  sputtering  features  followed  by  the  nuclei  and 
then  the  cytoplasm.  This  root  tip  section  was  0.5  ym  thick  and  took  —5  minutes 
to  completely  burn  through  to  the  substrate. 

Figure  6  presents  the  raw  <*°Ca+  ion  image  and  the  same  image  after  detector 
discrimination,  ion  intensity  conversion  and  differential  sputtering  corrections. 
Both  images  have  been  normalized  to  the  same  maximum  intensity.  Significant 
differences  in  the  distribution  of  calcium  in  the  tissue  are  observed  and  the 
corrected  image  more  accurately  indicates  the  presence  of  calcium  in  the  res¬ 
pective  features,  especially  in  the  cytoplasm.  This  image  more  closely  demon¬ 
strates  the  physiological  distribution  expected  in  the  in  vivo  situation. 

Quantitatively  there  has  been  a  change  in  the  relative  intensities  of  different 
features.  This  will  obviously  affect  the  results  when  converting  finally  from 
ion  intensity  to  concentration. 

The  relative  feature  sputtering  rates  R^  can  be  calculated  once  the 

burn-through  map  has  been  generated  by  ratioing  the  average  burn-through  rimes 

of  different  features.  Table  II  contains  the  average  R^  values  and  corresponding 

standard  deviations  for  cell  walls  and  nuclei  in  the  stele  cell  region  of 

radish  root  tip  cells  under  both  Ar+  and  02+  bombardment.  The  sputter  rates 
have  been  normalized  to  the  cytoplasm  sputter  rate,  the  slowest  sputtering 

feature. 

The  Rf  values  indicate-  that  differential  sputtering  makes  cell  walls 
appear  to  have  1.5  times  and  nuclei  1.2  times  more  calcium  (or  any  other 
element)  than  cytoplasm.  Like  results  have  been  observed  in  our  laboratory 
with  corn  and  pea'  root  tip  cells.  Galle  (9)  has  also  reported  nuclei  to 
sputter  faster  than  cytoplasm  in  frog  nucleated  red  blood  cells. 
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The  differentially  sputtering  nature  of  biological  tissue  may  be  due 
to  its  heterogeneous  chemical  composition.  In  the  case  of  plant  cells, 
cell  wall  is  primarily  composed  of  polysaccharides  (cellulose,  hemicell- 
ulose,  pectine,  etc.),  while  the  nucleus  is  basically  a  proteinous  material 
rich  in  DNA  (deoxyribonucleic  acid) ,  RNA  (ribonucleic  acid) ,  nucleoproteins , 
and  histone  proteins. 

A  bum-through  map  at  a  substrate  mass  also  revealed  that  epoxy 
(where  tissue  is  not  present)  was  the  last  to  sputter  through,  (i.e.,  after 
cytoplasm) .  This  helps  to  explain  the  influence  of  epoxy  in  the  differ¬ 
ential  sputtering  of  epoxy-infiltrated  biological  tissue.  In  the  case 
of  plants,  cell  walls  are  denser  material  than  nuclei  and  cytoplasm, 
hence  there  may  be  more  penetration  of  epoxy  (mass  per  unit  area)  into 

the  less  dense  cytoplasm  than  the  denser  nuclei  or  still  more  dense  cell  walls. 

It  is  interesting  to  note  that  both  primary  beam  species  induce 
similar  Rf  values  even  though  in  addition  to  the  differences  in  reactivity 
of  these  two  species  the  average  tissue  sputter  rate  for  an  Ar+  primary 
beam  is  greater  than  that  observed  with  02+  primary  beam. 

The  R^  values  indicate  the  possibility  that  our  assumption  of 
uniform  sputter  rate  with  depth  may  not  necessarily  be  valid,  at  least 
when  section  thickness  exceeds  0.5  (Jm.  The  decrease  of  the  Rf  value  with 
increasing  sample  thickness  indicates  significant  sample  homogenization. 

In  fact,  partial  diffusion  of  40Ca+  was  observed  to  occur  after  several 
minutes  of  bombardment,  possibly  due  to  local  sample  ^eating  or  charge 
induced  diffusion.  However,  the  basic  sample  structure  appears  unchanged 
as  cellular  features  are  well  defined  in  the  181Ta+  image  and  burn-through 
maps. 
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Factors  other  than  differential  sputtering,  such  as  sample  surface 
topology  (17)  and  matrix  effects  (18)  can  also  influence  the  observed 
secondary  ion  signal  intensity,  often  by  several  orders  of  magnitude.  Several 
methods  of  dealing  with  these  problems  have  been  examined  elsewhere  (19),  and 
are  beyond  the  scope  of  this  discussion.  While  the  influence  of  these 
factors,  particularly  the  matrix  effects,  may  be  larger  than  that  of  differential 
sputtering,  any  comprehensive  quantitation  procedure  will  have  to  include 
corrections  for  differential  sputtering. 

In  summary,  differential  sputtering  presents  a  problem  in  the  quanti¬ 
fication  of  ion  images,  particularly  for  biological  thin  sections  where  cellular 
organelles  differ  widely  in  their  chemical  composition.  Image  depth  profiling 
has  been  used  to  make  first-order  corrections  for  this  problem.  This  method 
of  correction  is  useful  when  a  thin  flat  sample  is  deposited  on  a  substrate. 

In  addition  to  biological  tissues  this  technique  should  also  be  applicable  to 
other  thin  films,  such  as  certain  types  of  semiconductor  devices  and  metall¬ 
urgical  samples. 
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Table  X.  Instrumental  Operating  Conditions 


Primary  Ions 
Primary  ion  energy 
Primary  ion  current 
Sputtered  area  of  sample 
Secondary  ion  polarity 
Transfer  optics 
Vacuus  at  sample 
Image  field-of-view 
Image  integration  period 
Time  between  images 


IQ  keV 

3.6x10  4  A/cm 
5 . 6x10  5  cm 
positive 
150  ym 
3x10  8torr 
50  ym 
2  seconds 
8-40  seconds 


Table  II.  Average  Rf  Values  by  Section  Thickness 


values 


Primary 

Beam 


Cellular 

Feature 


0.25  ym 
thick  section 


0.50  ym 
thick  section 


1.0  ym 

thick  section 


Rf  cell  walls 


1.51  ±  0.22 


1.46  ±  0.04 


1.19  ±  0.06 


Rf  nuclei 


1.20  ±  0.04 


1.16  i  0.02 


1.09  ±  0.06 


R^  cell  walls 


1.30  ±  0.08 


1.50  ±  0.02 


1.17  ±  0.01 


R_  nuclei 


1.18  ±  0.13 


1.18  ±  0.04 


1.15  ±  0.01 


I 


FIGURE  CAPTIONS 

Figure  1.  System  configuration 

Figure  2.  Light  micrograph  of  radish  root  tip  stele  cell  region 

(nu  »  nucleus ,  cw  -  cell  wall ,  cy  -  cytoplasm  ) 

Figure  3.  40Ca+  ion  micrograph  of  stele  cells  (nu  »  nucleus, 

cw  *  cell  wall,  cy  «.  cytoplasm) 

Figure  4.  181Ta+  ion  micrograph  of  substrate  150  seconds  into 

image  depth  profile 

Figure  5.  Time— domain  bum- through  map 

Figure  6.  Comparison  of  original  elemental  ion  image  and 
same  image  after  detector  discrimination,  ion 
intensity  conversion,  and  differential  sputtering 


corrections. 
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